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I. INTRODUCTION 

In practice, the control system designer rarely has the 
freedom to feed back the entire state of the system to 
achieve the desired system performance. Thus, an important 
problem from a practical standpoint is the determination of 
constant, output feed back gains for the control of systems 
with unmeasured states. 

The pole placement only and pole placement and robust- 
ness design procedure uses placement to establish a designer 
selected performance level and then a minimum singular value 
level to establish robustness. The pole placement only was 
relatively easy to implement through a numerical optimiza- 
tion routine. By using this numerical procedure it is also 
simple to incorporate robustness into the procedure along 
with the performance requirement. The technique, which 
utilizes a modern optimization routine, can significantly 
assist the designer in obtaining robustness in the face of 
cross coupling perturbations. It has been shown that the 
cross coupling problem can be detected by using classical 
open loop Bode diagrams as well as modern control analysis 
based on singular values. As currently employed, the pole 
placement and robustness program is used to obtain pole 
placement and robustness for a given set of starting gains 
and a selected optimization routine from the ADS program. 

The pole placement and robustness design procedure is a 
straight forward numerical optimization procedure for the 
practical application modern MIMO system analysis. The new 
aspects of the procedure are the implementation of both pole 
placement and robustness criteria within the same design 
program. 
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The pole placement and robustness design ro 
oped for this thesis has been used on X-22A A/C 
these studies the pole placement and robustness 
has proven capable of meeting the desired go 
placement and robustness and also brought to 
interesting aspects of the cross-coupling 
problem . 

The remainder of the thesis will present 
material on Robustness (SISO system) in Chap 
Robustness of MIMO systems in chapter Three. Po 
only design procedure in chapter Four, Pole p 
Robustness design in chapter Five, along with 
of Applications for X-22A a/c problem by inp 
function F*G analysis and output transfer functi 
ysis. Conclusions will be presented in the fi 
The computer programs used in the present an 
developed in [Ref. 9]. 
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II. ROBUSTNESS; SINGLE^INPUT S INGLE -OUTPUT SYSTEKS (SI SO) 



A review of the concept of robustness and 
the framework of a conventional SISO system w 
before pursuing the concepts in a more complic 
in the following chapters. A simple interp 
robustness is the ability of the system to tol 
perturbations. These perturbations could be in 
actuator failures, plant parameter uncertainty 
dynamics or nonlinear terms, or any one of 
perturbations to the nominal design of the syste 
The primary reason for feedback systems is 
of uncertainty within the system. By appropr 
feedback, properties that would lead to an uns 
may be controlled. When stability and robust 
are considered for a SISO system, frequency d 
concepts, using either Nyguist or Bode plots, 
used. Robustness in SISO systems is formulated 
the concept of gain and phase margins, both o 
readily available on the Nyguist or Bode diagram 
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1 





Figure 2.1 . Classical Bode Plot. 

Gain and phase margin can be defined in terms of the 
open-loop frequency domain plots in either the Bode or 
Nyquist format. Figure 2.1 depicts a classical Bode plot 
showing gain and phase margin determination from the plot. 
The Nyquist plot may also be used to obtain this informa- 
tion. Nyquist criterion states that if the open-loop 
transfer function G(s)H (s) has n poles in the right half 
plane and the limit of G (s) H (s) =constant as s — >oo then for a 
stable system the locus of G(s)H(s) will encircle the -1+jO 
point n times in the counterclockwise direction as s varies 
along the Nyquist contour. If there are no poles in the 
right half s plane then the locus will not encircle the -1 + 
jO point. The diagram in Figure 2.2 illustrates a nominally 
stable system. The gain and phase margin may be determined 
directly from the diagram.* 

Any change in the loop transfer function, provided the 
order of G(s)H(s) does not change, that changes the number 
of times the locus of G(s)H(s) encircles the (-1/0) point in 
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1 




Figure 2.2 Nyguist Plot of Stable System. 

the Nyguist plot causes the system to become unstable. This 
leads to the conclusion that the minimum distance of the 
locus of G(s)H(s) to the (-1/0) point is a measure of the 
system stability. This distance concept carries over 
directly to the Multi-Input Multi-Output (MIMD) system as 
will be shown in the next chapter. Examples of a multipli- 
cative perturbation and an additive perturbation illustrate 
this idea. Figure 2.3 is an additively perturbed system. 
Figure 2.4 shows the Nyguist plot for this system. Assuming 
that the plant is itself stable and the perturbations are 
also stable the diagram may then be used to determine how 
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Figure 2.3 Additively Perturbed System. 

near the system is to instability for the given perturba- 
tion . 

Since the system is stable the (-1,0) point is encircled 
the correct number of times by the nominal plant. If the 
locus of g(jOU) in the diagram is warped until it passes 
beyond the (-1,0) point then clearly the number of encircle- 
ments of this point will change and the system will become 
unstable, assuming the order of the plant is not changed by 
the perturbation. To keep the locus of points from moving 
beyond the (-1,0) point equation 2.1 must hold. 
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Figure 2.4 Additive Nyquist Plot. 



1 Ag(jto) l<I 1 + g (ju)) I (2.1) 

This condition is illustrated in Figure 2.5. The right-hand 
side of equation 2. 1 is just the magnitude of the return 
difference transfer function of the nominal system. The 
multiplicative case is depicted in Figure 2.6 with its asso- 
ciated Nyquist plot in Figure 2.7. The requirement for 
stability is similar to the additive case and may be stated 
in equation 2.2 

1 Ag(jw) l < I i ♦ (g(jco) -*) i 



1 6 



( 2 . 2 ) 



1 



Im (Jto) 




Figure 2.5 Hyguist for Inequality Additive Condition. 



The above arguments will be applied again in Chapter 3 to 
develop multivariable stability and robustness properties. 
With this basic review of the concepts of stability and 
robustness in the classical SISO system complete, the next 
chapter will extend some of these basic concepts to the dIMO 
system. 
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III. HULTIVABIABLE SYSTEMS ROBUSTNESS 



A generalization of the SISG Hyquist theory discussed in 
the previous chapter has been made for the MIMO problem. 
This generalization leads directly to the apt- licatioi. of 
singular value concept. The generalization is expressed in 
the form of the multivariable Nyguist theorem which requires 
that a closed loop stable system have the same number of 
counterclockwise encirclements of the origin by the locus of 
the det (I+G (jio) ) as the number of open loop poles that are 
in the right half plane. This theorem is formally stated as; 
let N[f (s) ] denote the number of clockwise encirclements of 
(-1,0) by the locus of £ (s) as s traverses the contour D 
of Figure 3.1 in a clockwise sense. The closed-loop system 
will be stable if and only if for all R sufficiently large 
N[f (s) ] = -P 

where N = number of encirclements 

f (s) = -1 + det[ I+G (s) ] = (s) / <f ct _ (s) -1 and 

P = the number of closed right-half plane zeros of 

9 ., < s > • 

The application of the Nyguist theorem comes through the 
fact that a multivariable system will not be robust to 
modelling errors if the return difference matrix, I+G, is 
nearly singular for some frequency. If I + G is nearly 
singular a small change in G may make I + G exactly 
singular. This causes the det (I + G) to become zero and the 
Nyguist encirclement count to change indicating an unstable 
system. It is possible for very small changes in I+G to 
produce large changes in the determinan tof 1+3. The matrix 
I+G 
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Figure 3.1 Nyguist D Contour. 

/ 

1 0/s+a 9.9 9/s«- a ' 

10/s+a 10/s+a 

has determinant 0.1/ (s+a) 2 . If the element p is changed 
by only one percent to 9.9/s+a the determinant becomes 
1.1/ (s+a) 2 which is a significant change in the determinant- 
value. Therefore, it is evident that det (I + G) is not an 
accurate measure of how near the return difference is to 
singularity. Researchers, in the field of controls 
[Ref. 1], [fief. 2], [Bef. 3], [Ref. 4] have used singular 
value analysis to determine how near the return difference 
matrix is to singularity. 

Since the number of encirclements of the Nyguist diagram 
changes as f (s) passes through the -1 point or when det(I+G) 
is zero it is important to find how near the return differ- 
ence matrix I+G is to being singular. This nearness to 
singularity can be interpreted as the distance of matrix I+G 
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to the critical point, -1. A guantity which can be used to 
express the nearness to singularity of the matrix is the 
minimum matrix singular value denoted by^ . Given a matrix 
A the singular value may be expressed by eguation 3. 1 

Q\ (A) = min ( A* ( A* A) ) (3.1) 

A 

u 

where ^ (A_ A) is the eigenvalue of the complex conjugate 
transpose of A_ times A. A basic MIMO linear system is 










Figure 3.2 Basic Multi-input Multi-output System. 

depicted in Figure 3.2. An additive perturbation to the 
plant is shown in Figure 3.3. If the plant is stable before 
the perturbation is added to the system the Nyguist theorem 
will be satisfied and the locus of GH will not encircle the 
-1,0 critical point. When the perturbation is added to the 
system as long as the Nyguist locus is not forced to 
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encircle the -1 point the system, will remain stable. A 
sufficient condition, recalling the SISO discussion in 
chapter 2, for the perturbed Nyguist plot not to change 
encirclements is that the norm of the perturbation Ag 
remain less than the norm of the return difference matrix as 
expressed in equation 3.2. 

II A GJjWJ) I ! < VI I (I+G) -i) I I (3.2) 

The condition u> £0 will that the locus of the det (I+G) 
does not pass through the -1 point. If the 1 or Euclidean 
norm is assumed for this condition the equation 3.2 may be 
expressed in terms of singular values as equation 3. 3. 
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CN (AG) < Qs (I + G) 



(3.3) 



This result states that as long as the maximum singular 
value of the pert ur bation matrix A G is below the minimum 
norm value of the return difference matrix the system will 
remain stable. The problem of guaranteeing robustness 
becomes that of finding the largest norm of the perturbation 
guantity, the largest singular value, for which the smallest 
norm or singular value of the return difference matrix will 
remain non-singular. 

The multiplicative form for a system such as Figure 3.4 
gives the similar norm equation in equation 3.4 

1 1 AG (j£U) 1 1 <1/1 I (I + (G)-»)-MI (3.4) 

03 >0. which may be expressed as 

(N ( AG) < 0\ ( I + G-» ) (3.5) 

Singular value decomposition software is readily available 
to determine how near the matrix I+_G or 1+ (G) -1 is to 

singularity. 
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IV. POLE PLACEMENT ONLY DESIGN PEOCEDDBE 

The designer may use either objective, constraint or a 
combination of functions to secure the desired pole loca- 
tions. As currently implemented in the program the cost or 
objective portion of the pole placement procedure is 
constructed as equation 4.1 

i 

OBJ = YL ( Am- A r . ) * + ( A,..- Xi, ) 2 < t *- 

A =1 ' 

vhere A R = real eigenvalue 

Ax = imaginary eigenvalue 
A*p = desired' eigenva lue location 
7\jd = desired eigenvalue location 
The constraint formulation is a function that must be 
kept negative or the constraint is violated. It is written 
as equation 4.2 

g(j) = \j (Awr 2 + ( Ajd" A r; ) 2 - r (4.2) 

where r is a tolerance circle established as a function of 
pole placement position. Since the aim of the optimizer is 
to keep g negative any time the function of the 

constraint is greater than r the constraint will become 
active, i.e. violated. The optimizer will then attempt to 
move the constraint to the inactive status by adjusting the 
design parameters of the system. For obtaining the feedback 
gains and desired poles use the 'CONXSV* program [Ref. 9], 
weight function 1 equal 1 the other weight functions are set 
to 0. After obtaining the feedback gains the NPS ’OPTSYS* 
program was used to plot the Bode, Nyquist diagram and pole 
zeros plot. 
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V. POLE PLACEMENT AHD ROBUSTNESS DESIGN PROCEDURE 



Consideration of implementation of the frequency domain 
or robustness portion of the design procedure begins with 
the concept of tllMO phase and gain margins. Several useful 
theorems on singular value analysis of multiloop systems are 
presented in [Ref. 3]. One of these theorems relates the 
matrix singular value of the return difference function to a 
para meter, qC , and further shows that as long as the maximum 
singular value of the perturbation function (L _1 - I) 

remains less than thisoC, the system reaains stable. The 
value of oi is then related to gain and phase margins of the 
MIMO system. The relationship developed is given in equa- 
tions 5 . 1 and 5. 2: 

gain margin = GM =1/(1 +d Q ) (5.1) 

phase margin = PM = +cos -1 ( 1 - d t 2 /2) (5.2) 

provided that equation 5.3 holds. 

(I + G) > oL 0 
for some c(p S 1 

These phase and gain margins are guaranteed in 
simultaneously . 

Universal gain and phase margin curves, [Ref 
on the minimum singular values of the return 
matrix are developed from equation 5.4. 

(L-i-I) =max (1- 1/k*) 2 + 2/k n (1-sosft) (5.4) 



(5.3) 

every loop 

. 5], based 
difference 
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foe all n with k^ > 0. These curves shown in Figure 5. 1 
allow the designer to pick a singular value that corresponds 
to a specific gain and phase margin for a given system. 




Figure 5.1 Universal Gain. and Phase Singular Value Plot. 

Since the universal curve in Figure 5.1 provides a 
convenient method of specifying gain and phase margins in 
terms of singular values the robustness portion of the pole 
placement and robustness design procedure uses the minimum 
singular value level of the return difference matrix to 
determine the robustness. The minimum singular value level 
is formulated as an objective or constraint function in 
equation 5.5 

J= E(»ax (0, ((£- <r (jui,p)) ) 2 (5.5) 

4 
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The optimization procedure may be used to change feedback 
gains until the minimum singular value is raised above this 
desired design level. Although the same formulation can be 
used as a negative constraint function it has not been 
implemented as such within this program. There are numerous 
ways the singular value formulation could be implemented 
within the program by changes of the code if design require- 
ments forced such changes. 

The pole placement and robustness design program is 
based on the ADS code to implement the design variable 
selection procedures. The pole placement and robustness 
program is used to provide designs for state or output feed- 
back problem. 

For the state or output feedback design program the user 
must input the plant matrices A.,B, C_ and initial starting 
values for the feedback matrix F. The matrices correspond to 
the following linear differential system: 



X = Ax_ + 


Bu 


(5.6) 


y = Cx 




(5.7) 


u = -Fx 




(5.8) 


As the 


design program is currently coded 


the user may 


run output 


feedback or state feedback design 


by specifying 



the C matrix as the diagonal(I) matrix for state feedback. 
The program relies on initial starting values of the feed- 
back gains, F. 

The ability to select acceptable starting values for the 
feedback gains will make the procedure more efficient in 
operation. As currently employed, the program is used to 
obtain pole placement and robustness for a given set of 
starting gains and a selected optimization routine from the 
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ADS program. If the optimizer is not able to meet the 
desired design goals on this program run two options are 
available. First, change to a different optimization routine 
from the list of available ADS routines and rerun the 
problem. This was usually successful in improving the 
design. Second, the designer uses a new set of starting 
values for the feedback gains and repeats the design proce- 
dure. Both options might be used on particularly difficult 
cases . 

The pole placement and robustness design algorithm 
computes input additive, output additive, input multiplica- 
tive, and output multiplicative singular values. The versa- 
tility of the pole placement and robustness design is 
obtained by incorporating a state of the art optimizer 
routine ADS, with currently available singular value compu- 
tation routines. 



30 



VI. APPLICATIONS (X22 A LONGITODIN AL PROBLEM) 



This chapter will deal with a more practical application 
of the numerical optimization program. In this problem the 
combined pole placement only/ robustness design procedure 
will be applied to the linear longitudinal dynamic channels 
of an X22A V/STOL a/c [Ref. 8]. The model is that of the 
longitudinal dynamics of an X22A V/SIOL a/c at low altitude 
and airspeed = 65 knots. The dynamic model of tne system is 



x = A x + B u 



( 6 . 1 ) 



X ={u,w r g,0) 



( 6 . 2 ) 



a = (<£<<£> 



(6.3) 



where 



A = 



and 



' -0.18 -0.03 

-0.2 -0.55 

- 0.01 - 0.02 
0 0 
where 

' -0.36 0.52 



B 



0 . - 1 . 



9.57 -31.87 

109. 43 2. 78 

- 0 . 1 0 

1 . 0 



0.33 0.02 

0 . 0 . 

^ J 

with full state available for feedback. Table 1 is a summary 
of parameters. The system is not open-loop stable. 

The control law is formulated to satisfy the desired 
performance specif ication. Equation 6.4 is the basic control 

law . 



TABLE 1 



X22A V/STOL a/c Parameter Definitions 



Variable 


Units 


Description 


0 


f t/sec 


Body axis forward velocity perturbati) n 


W 


f t/sec 


Vertical velocity perturbation 


Q 


rad/sec 


Pitch rate 


0 


rad 


Pitch attitude 


de 


inch 


Elevator position 




inch 


Throttle position 



u = -Fx 



(6.4) 



The design studies presented up to this point have been 
based on breaking the system locp at the input as shown in 




Figure 6.1 System Block Diagram. 



32 



Figure 6.1. In multivariable tneory the location of the 
break in the loop changes the rsturn difference for the 
system and the transfer function formulation. In Figure 6.1, 
number 1 depicts a system with an input loop break point 
while number 2 depicts an output loop break point for output 
return difference determination. The return difference func- 
tion for point 1 is written as I+F3 while the return differ- 
ence for point 2 is I+GF. In figure number 1 the transfer 



function FG 


is 




FG = F’ * (SI - A) i *E 




F’ = F*C 


and GF is 


GF = C* (SI - A) 1 *B ’ 




B» = B*F 



A. INPUT TRAHSFEB FUNCTION F*G ANALYSIS 

Using the NPGS ’CONXSV’ program [Ref. 9], to get the 
desired feedback constants, the pole placement only program 
was first run to get the desired pole and the computed 



TABLE 2 



Pole Only Desired Pole and Computed Pole 



Desired pole 


Computed pole 


-1.2 1 . 6 j 
-1.2 -1.63 
-1.0 J 
-0.5 


-1.19889 1.59415i 

-1.19889 -1.594151 

-0. 98902 J 

-0. 49770 



pole (see Table 2) . The feedback matrix F is 



F 



-0. 10853 0.50675 13. 97235 

0.19809 -0.10997 -90.41902 



7.73604 

12.42874 



33 




Figure 6.2 Pole Placement Only Singular Value Plots. 



the objective function • is 0.0002 and the additive minimum 

N. 

singular value vs frequency is seen in Figure 6.2. This 
plot shows that pole placement only design has very low 
minimum singular values for the return difference matrix. 
Pole placement only goes as low as -22 db or 0.0827 rad near 
2 rad/sec in frequency. Using the universal gain and phase 
diagram as discussed in Chapter 5 this equates to a gain 
margin of about 0.95 db to 1.15 db and a phase margin of 4.7 
degrees. These phase and gain margins are quite small, 
showing the need to run a pole placement and robustness 
design program. It is assumed that the eigenvalues (pole 
locations) of the system, as developed in [Eef. 7], are the 
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required poles for the performance criteria. Once the pole 
locations are set, robustness criteria must be selected. 
From the universal gain and phase margin curve discussed 
earlier, singular value levels were selected for this 
problem. The singular value level chosen was 0.6 rad. This 
corresponds to a gain margin of -4.0 db to 8 db and a phase 
margin of about 35 degrees. 

For a singular value level cf 0.6 rad the pole placement 
and robustness design routine places the poles as shown in 
Table 3. The slight differences in these pole locations 



TABLE 3 



Pole Placement and Robustness Design Pole Locations 



Desired pole 


Computed pole 


- 1.2 1 . 6 j 

- 1.2 -I .63 

- 1.0 
-0.5 


-1. 19297 1 . 59533 j 

-1. 19297 -1. 595333 

-0. 99886 
-0. 45460 



have an insignificant effect on the performance. 
The feedback constant F is 



F 



-0.07493 -0.08746 8.0529 16.30763 

0.27085 -0.23019 -2.2919 37.66504 



The objective function is 0.0000 and the additive 
minimum singular value vs frequency is seen in Figure 6.3. 
This plot shows that in the robustness design the gain 
adjustment moves the minimum singular value from about 0.08 
rad, with very poor phase and gain margins, to a level of 
0.645 rad, above the desired values of gain and phase. 
After obtaining the feedback gains the NPGS 'OPTSYS' program 
was used to obtain the necessary data, time response plot, 
Eode plot, and pole-zero map. The pole placement only 
design closed loop time response plot is shown in Figure 
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Figure 6.3 Robustness Design Singular Value Plots. 

6.4. It shows good pitch rate and pitch attitude responses 
to <£ step input, but slow responses to the forward 
velocity perturbation and vertical velocity perturbation to 
eft step input. The input transfer function F*G robustness 
design time response plot is shown in Figure 6.5. This 
indicates a better time response for pitch rate and pitch 
attitude inputs than pole placement only design, also good 
vertical velocity perturbation response to gQ step input, 
but sluggish response for the forward velocity perturbation 
[Ref. 8]. Comparing the loop Eode plots of the pole place- 
ment and robustness design transfer function, F*G, shows an 
increase in robustness. 
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Step Input Transient Response 




q/JL Step Input Transient Response 




Figure 6.4 Pole Placement Only Design Time Response. 
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Step Input Transient Response 




of Step Input Transient Response 
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Figure 6.5 F*G Robustness Design Time Response 
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Figure 6.6 Pole Placement Only Design F*G 1:1. 
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Figure 6-7 Robustness Design F*G 1:1. 
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Figure 6.8 Pole Placement Only Design F*G 1:2. 
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Figure 6-9 



Robustness Design 



F*G 1: 2. 
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Figure 6.10 Pole Placement Only Design F*G 2:1. 
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Figure 6.11 Robustness Design F*G 2:1. 
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Figure 6.12 Pole Placement Only Design F*G 2:2. 
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Figure 6.13 Robustness Design F*G 2:2. 
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FREQUENCY (w, rad/sec) 



The F*G transfer functions bode plot numerical results are 
shown in Table 4. 



TABLE 4 



F*G Transfer Functions Bode Plot Besults 



Channel 


Bandwidth 
rad/sec 
pole SV 


PM 

deg re e 
pole SV 


GM 

Db 

pole 


SV 


Gain 

Db 

pole 


SV 


1 


1 


6. 1 


3.2 


58 


80 


• 


• 


34 


25.7 


1 


2 


30 


2.83 


89.7 


-2 0 


-32 


-7.2 


28 


27.6 


2 


1 


0. 3 


0.19 


9.5 


53 


8. 8 


• 


19.4 


8.6 


2 


2 


0. 17 


0 .64 


170 


106 


- 7.3 


• 


13. 5 


14. 1 



The most noticed change in the overall system, however, 
is in the transfer function input 1 output 2 . This is the 
change at the minimum singular value position, which greatly 
reduces the gain and bandwidth in this channel through a 
change in feedback gain. The optimizer routine nrings the 
entire system gains to more balanced conditions and recovers 
a highly robust design. 

The gain changes' associated with the robustness improve- 
ment cause the zeros of the various closed-loop pole-zero 
diagram of the closed-loop transfer matrices to move. A 
comparison of the eight pole-zero diagrams is shown in 
Figures 6.14 to 6.17. 

The significant feature of these pole-zero diagrams is 
the shift in the zeros of the optimized design in a direc- 
tion that attempts to equalize or balance the frequency 
response for frequencies in the vicinity of the minimum 
singular values. The pole-zero diagram of channel input 1 
output 2 will be discussed as an example of this effect. In 
figure 6.12 the pole only design zeros are located at about 
-0.31 and -0.78. When the pole placement and robustness 
routine has completed the feedback gain modification these 
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Figure 6.14 Pole-Zero Map for 1:1 and 1:2 
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Figure 6.15 Pole-Zero Map for 1:3 and 1:4. 
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Figure 6.16 Pole-Zero Map for 2:1 and 2:2. 
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Figure 6.17 Pole-Zero Hap for 2:3 and 2:4. 
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zeros have shifted to -0.17 and -0.56. The effect of these 
zero shifts is to combine with the pole, locations to 
equalize the frequency response and increase DC gain of 14 
db for the same bandwidth, at the minimum singular value 

freguency 1 rad/sec as depicted in Figures 6.18 and 6.19. 

Zero shifts for the remainder of the transfer functions 
provide similar results in the other channels. By moving 
toward the frequencies associated with the minimum singular 
values the zeros have balanced the overall freguency 

response of the system in each channel. While the channel 
gain modification is the primary mechanism for robustness 
recovery, the zero shift associated with the feedback gain 
changes is directly related to the overall freguency 

response of the system. A robustness design meeting the 
required pole locations and a robustness singular value 
level of 0.6 rad provides adeguate gain and phase margin for 
the design. 



B. 



ODTPDT TRANSFER FUNCTION G*F ANALYSIS 



Using the 
desired output 



NPGS 1 CO H X S V 1 program [Ref. 9], to get the 
transfer function feedback constants, the 
pole placement only program was first run to get the desired 
pole and the computed pole (; 



matrix F is 
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Figure 6.18 Pole Placement Only Design F*G Closedloop 1:2. 
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Figure 6.19 Robustness Design F*G Closedloop 1:2. 
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FREQUENCY (cj, rad/sec) 





Figure 6,20 Pole placement only G*F Singular Value Plots. 

matrix. Pole placement only goes as low as -48 db or 0.004 
rad near 1 rad/sec in frequency. Using the universal gain 
and phase diagram as discussed in Chapter 5 this equates to 
a gain margin of about -0.003 db to 0.003 db and a phase 
margin of 0.5 degrees. These phase and gain margins are 
quite small, showing the need to run a pole placement and 
robustness design program. It is assumed that the eigenva- 
lues (pole locations) of the system, as developed in 
[Ref. 7], are the required poles for the performance 
criteria. Once the pole locations are set (see Table 5), 
robustness criteria curve discussed earlier, singular value 
levels are made for this problem. The singular value level 
chosen is 0.6 rad. This corresponds to a gain margin of -4.0 
ab to 8 db and a phase margin of about 35 degrees. 
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For a singular value level of 0.6 rad the pole placement 
and robustness design routine places the poles as shown in 
Table 5. There is a trade off between performance and 
robustness in many control problems. IN output transfer 
function G*F singular value calculation# the shift in this 
pole can have a significant effect on the performance, it 
was found necessary to move the pole at -0.5 to -0.05 in 
order to achieve the desired robustness. The slight differ- 
ences in these pole locations an insignif icant effect on the 



TABLE 5 

G*F Pole Placement and Robustness Design Pole Locations 



Desired pole 


Computed pole 


-1.2 1.61 
-1.2 - 1 . 6 j 

-1.0 
-0.05 


-1.20944 1 . 58904 j 

-1.20944 -1.589041 

-1 . 02955 
-0. 05336 



performance. The feedback matrix F is 

F = f -0.27091 0.18869 0.0 0.0 

1.88298 -1.60017 0.0 0.00914 

V. 

The objective function is 0.0002 and the 



additive 



output singular value vs frequency is seen in Figure 6.21. 
This plot shows that in the robustness design the gain 
adjustment moves the minimum singular value from about 0.004 
rad, with very poor phase and gain margins, to a level of 
0.612 rad, close enough to the desired values of gain and 
phase. After obtaining the feedback gains the NPGS ’CPTSYS 1 
program was used to obtain the necessary data. Bode plot, 
pole-zero map, and time response plot. The G*F robustness 
design time response plot is shown in Figure 6.22. The 
response of pitch rate is relatively good but sluggish 
response of pitch attitude to c£ step input. The performance 
of the pitch attitude was degraded. It is possible that a 
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Figure 6.21 Eobustness Design G*F Singular Value Plots. 

further change in the pole location could achieve better 
performance and still meet the robustness criteria. The 
body axis forward velocity and vertical velocity perturba- 
tion responses are the same as in the pole placement only 
design(see Figure 6.4). Comparing the loop Bode plots of 
the pole placement and robustness design transfer function, 
G*F, shows an increase in robustness. 

Considering a SISO Bode analysis, the channel input 1 
output 1 bode plot are shown in Figure 6.23 and 6.24. The 
bandwidth (BW) shows a small increase, 0.6 rad/sec, and the 
phase margin decrease of 26 degrees. The DC gain drops from 
21.4 db to 17.3 db, a 4 db decrease and the slope change is 
similar to the pole only design as seen in Figures 6.23 and 
6.24. This channel only shows a smooth curve. 
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Figure 6.22 G*F Robustness Design Time Response. 
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Figure 6.23 Pole Placement Only Design G*F 1:1. 
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Figure 6.24 Robustness Design G*F 1:1. 
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FREQUENCY (cj, rad/sec) 



Considering the channel input 1 output 2 with bode plots 
shown in Figures 6.25 and 6.26, a shift in the bandwidth (3W) 
from 2.4 rad/sec to 3.2 rad/sec is shown. The DC gain 
increases from 12.7 db to 15.8 db at the minimum singular 
frequency of 1 rad/sec and an increase in the phase margin 
of 28 degrees. Figures 6.25 and 6.26 show the similar 
frequency response curve, slope changes by zeros move the 
one pole position 0.8 rad/sec for pole-zero cancel and the 
other one zero move to close the minimum singular value 
frequency 1 rad/sec and two pole location. The increase in 
the DC gain and phase margin yields an increased tolerance 
to perturbation. The channel input 2 output 1 bode plots 
are shown in Figures 6.27 and 6.28. The DC gain decreases 
16 db and the bandwidth shifts left 0.8 rad/sec, increasing 
the phase margin 151 degrees. The zero shift has the effect 
of smoothing the frequency response curve in the vicinity of 
the frequency of the minimum singular value. The channel 
input 2 output 2 bode plots are shown in Figures 6.29 and 
6.30. These figures show a small shift left in the band- 
width from 4.5 rad/sec to 2.8 rad/sec. The DC gain decreases 
12 db with an increased phase margin of 56 degrees. This 
results in a smooth curve, caused by the zero location shift 
in the same manner as previous channels (channels 
1; 1, 2; 1 , 2; 2) . The G*F transfer functions bode plot numer- 
ical results are provided in Table 6. 

The noticeable change in the overall system, however, is 
in the transfer function input 1 output 2 . This is the 
change at the minimum singular value position, increasing 
the gain and phase margin through a change in feedback gain. 
The optimizer routine brings the system gains to a more 
balanced condition and recovers a highly robust design. 

The gain changes associated with the robustness improve- 
ment cause the zeros of the various closed-loop pole-zero 
diagram of the closed-loop transfer matrix to move. A 
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Figure 6.25 Pole Placement Only Design G*F 1:2. 
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Figure 6.26 Robustness Design G*F 1:2. 
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Figure 6.27 Pole Placement Only Design G*F 2:1. 
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Figure 6.28 Robustness Design G*F 2:1. 



66 



fact tvsova 



(38P) 3SVHd 




Figure 6.29 Pole Placement Only Design G*F 2:2. 
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Figure 6.30 Robustness Design G*F 2:2. 



68 



FREQUENCY (a>, rad/sec) 



TABLE 6 



G*F Transfer Functions Bode Plot Results 



channel 


Bandwidth 
rad/sec 
pole SV 


PM 

degree 
pole SV 


GM 

Db 

pole 


SV 


Gain 

Db 

pole 


SV 


1 ; 1 


0.95 


1.56 


49 


23 


3.9 


• 


21.4 


17.3 


1 : 2 


2.36 


3.2 


211.3 


239 


• 


• 


12.6 


15.8 


2 : 1 


2.2 


1 .43 


80 


231 


m 


• 


37.3 


21.5 


2 : 2 


4. 5 


2.8 


10.7 


17.5 


• 


• 


29.1 


17.5 



comparison of the eight pole-zero diagrams is shown in 
Figures 6.31 to 6.34. 

The significant feature of these pole-zero diagrams is 
the shift of the zeros of the optimized design in a direc- 
tion that attempts to. equalize or balance the frequency 
response for frequencies in the vicinity of the minimum 
singular values. The pole-zero diagram of channel input 1 
output 2 will be discussed as an example of this effect. In 
figure 6.30, the pcle only design zeros are located about 
-0.31 and -0.78. When the pole placement and robustness 
routine completes the feedback gain modification these zeros 
shift close to -1.7 and 0.28. The effect of these zeros 
shifting is to combine with the pole locations to equalize 
the frequency response and increase the DC gain 25 db at the 
same bandwidth and the minimum singular value frequency 1 
rad/sec, as depicted in Figures 6.35 and 6.36. Zero shifts 
for the remainder of the transfer functions provide similar 
results in the other channels. By moving toward the 
frequencies associated with the minimum singular values the 
zeros have balanced the overall frequency response of the 
system in each channel. While the channel gain modification 
is the primary mechanism for robustness recovery the zero 
shift associated with the feedback gain changes is directly 
related to the overall frequency response of the system. 
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Figure 6.31 Pole-Zero Hap for 1:1 and 1:2. 
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Figure 6.32 Pole-Zero Map for 1:3 and 1:4. 
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Figure 6.33 Pole-Zero Map for 2:1 and 2:2. 
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Figure 6.34 Pole-Zero Map for 2:3 and 2:4. 
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Figure 6.35 Pole Placement Only Design G*F Closedloop 1:2. 
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Figure 6.36 Robustness Design G*F Closedloop 1:2. 
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Figure 6.37 Robust Design F*G Output Singular Value Plots. 
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VII. CONCLUSIONS 

The pole placement and robustness design routine coupled 
with the Automated Design Synthesis program provides the 
designer an excellent tool with which to attack the robust 
design problem. The pole placement and robustness design 
routine has demonstrated the capability of providing designs 
that solve the X22A V/STOL A/C longitudinal dynamic problem 
caused by cross-coupling perturbations which reduce robust- 
ness in multivariable systems. This design improvement is 
accomplished by modifying the system feedback gains in such 
a manner that the gain in channels that are affected by 
cross-coupling perturbations is equalized with other system 
gains to reduce this cross-coupling effect. The gain changes 
are accompanied by zero shifts which also influence the gain 
distribution and freguency response of the system. 

The lack of robustness can be discovered in two ways. 
The first method is to plot the open-loop Bode plots of each 
element of the transfer matrix and look for extremely high 
gains or bandwidths relative to the other transfer func- 
tions.' The second method examines the singular values of 
the return difference matrix for magnitude. The joint anal- 
ysis method tells where the robustness problem occurs in the 
design. Low singular values correspond to low robustness. 
The pole placement and robustness design routine can 
increase robustness by modifying feedback gains to reduce 
the effect of cross-coupling within the system. Observing 
the gain modification made by the pole placement and robust- 
ness routine the critical channel within the system that 
affects the robustness can be determined from the 3ode 
plots. The pole placement and robustness routine feedback 
gain changes also cause zero shifts during the robustness 
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recovery. The gain on the open loop Bode plot for the 
affected cross-coupling channel is adjusted and the closed- 
loop zeros, as seen on the pole-zero diagram, are shifted. 
This zero shift is in a direction which will combine with 
system poles to smooth the frequency response diagram in the 
vicinity of the minimum singular value. The performance of 
the robustness design input transfer function F*G provides 
the best response, but the output transfer function G*F 
indicates the performance of the pitch attitude is degraded. 
It is possible that a further change in the pole location 
would provide better performance and still meet the robust- 
ness criteria. 
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